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540 FORM VOLUMETRIC ELEMENT GROUPS 

541 SELECT ONE VOLUMETRIC ELEMENT PER GROUP 

522 DETERMINE LOCAL TEMPERATURE FOR SELECTED VOLUMETRIC 
ELEMENTS, BASED ON H0P1H1 (T)+P2H2 (T)+P3H3 (T) 

523 DETERMINE AT THAT TEMPERATURE GIBBS FREE ENTHALPIES OF PHASES 
AND LEVELS OF CONCENTRATION OF PHASES AT THE PHASE BOUNDARIES 

524 DETERMINE THE NUMBER OF PHASES TO BE OBSERVED 
S RESOLVE THE DIFFUSION EQUATION 

528 AUSTEN ITE? 

529 RESOLVE STEFAN PROBLEM FOR AUSTEN ITE-FERRITE AND OTHER PHASE 

530 RESOLVE STEFAN PROBLEMS FOR AUSTENITE-FERRITE AND FOR 
AUSTENITE-CEMENTITE: OBSERVE STEFAN PROBLEM COUPLING 

537 DETERMINE NEW BOUNDARY PHASE POSITION BASED ON AX 

538 P3(ZF) = NON LINEAR FUNCTION (POSITION) 

542 RESTART RESULT FOR OTHER VOLUMETRIC ELEMENTS 

543 RESOLVE THERMAL CONDUCTION EQUATION FOR ALL VOLUMETRIC 

ELEMENTS 



(57) Abstract: The present invention concerns a computer (4) which determines a 
subsequent state (ZF) of a steel volume (1), based on an instantaneous initial state 
(ZA) of said steel volume (1) and at least one volumetric surface, the temporary in- 
fluence quantities (W) acting on said steel volume (1), by resolution of an equation 
of thermal conduction and phase change. The states (ZA; ZF) include for at least 
one volumetric element (9) of the steel volume (1), a local distribution in concen- 
tration (K) of a alloy element mobile in the steel, the local proportions (pi, p2, p3) 
of the modelled phases of the steel and a quantity (H) describing the local energy 
content of the steel. The phases comprise austenite and another phase, generally, 
ferrite or cementite. In the context of the change equation, the concentration levels 
(kl, k3; k2, k4) of the mobile alloy element, which are located on either side of the 
phase boundary (11, 12), between the austenite and the other phase are determined. 
Moreover, it is determined, by resolution of the Stephan problem, whether and how 
the distribution in concentration (K) of the mobile alloy element is modified in the 
austenitic zone of the volumetric element (9) observed and whether and to what 
extent (8x\ 8x") the phase boundary (11, 12) is thereby displaced. The proportions 
(pi, p2, p3) of phases are then determined based on the position of the limite phase 
(11, 12). 

(57) Zusammenfassung: Ein Rechner (4) ermittelt anhand eines momentanen 
Anfangszustands (ZA) eines Stahlvolumens (1) und mindestens einer uber 
die Volumenoberflache auf das Stahlvolumen (1) einwirkenden momentanen 
Einflussgrosse (W) durch Losen einer Warmeleitungsund einer Phasenumwand- 
lungsgleichung einen Folgezustand (ZF) des Stahlvolumens (1). Die Zustande 
(ZA, ZF) umfassen fiir mindestens ein Volumenelement (9) des Stahlvolumens 
(1) eine lokale Konzentrationsverteilung (K) eines im Stahl beweglichen 
Legierungselements, lokale Anteile (pi, p2, p3) von modellierten Phasen des 
Stahls und eine den lokalen Energieinhalt des Stahls beschreibende Grosse 
(H). Die Phasen umfassen Austenit und eine weitere Phase, in der Regel Ferrit 
oder Zementit. Im Rahmen der Umwandlungsgleichung wird ermittelt, welche 
Konzentrationen (kl, k3; k2, k4) des beweglichen Legierungselements beidseits 
einer Phasengrenze (11, 12) 

[Forts etzung auf der ndchsten Seite] 



WO 2005/076092 Al I II III II II 11 Hill 1 1 1 II 1 1 III II 1 1 II I II 1 1 1 1 III ill I III lllllll lllll llll 



(81) Bestimmungsstaaten (soweit nicht anders angegeben, fur 
jede verfiigbare nationale Schutzrechtsart): AE, AG, AL, 
AM, AT, AU, AZ, BA, BB, BG, BR, BW, BY, BZ, CA, CH, 
CN, CO, CR, CU, CZ, DE, DK, DM, DZ, EC, EE, EG, ES, 
FI, GB, GD, GE, GH, GM, HR, HU, ID, IL, IN, IS, JP, KE, 
KG, KP, KR, KZ, LC, LK, LR, LS, LT, LU, LV, MA, MD, 
MG, MK, MN, MW, MX, MZ, NA, NI, NO, NZ, OM, PG, 
PH, PL, PT, RO, RU, SC, SD, SE, SG, SK, SL, SY, TJ, TM, 
TN, TR, TT, TZ, UA, UG, US, UZ, VC, VN, YU, ZA, ZM, 
ZW. 

(84) Bestimmungsstaaten (soweit nicht anders angegeben, fur 
jede verfiigbare regionale Schutzrechtsart)'. ARIPO (BW, 
GH, GM, KE, LS, MW, MZ, NA, SD, SL, SZ, TZ, UG, 



ZM, ZW), eurasisches (AM, AZ, BY, KG, KZ, MD, RU, 
TJ, TM), europaisches (AT, BE, BG, CH, CY, CZ, DE, DK, 
EE, ES, FI, FR, GB, GR, HU, IE, IS, IT, LT, LU, MC, NL, 
PL, PT, RO, SE, SI, SK, TR), OAPI (BF, BJ, CF, CG, CI, 
CM, GA, GN, GQ, GW, ML, MR, NE, SN, TD, TG). 

Veroffentlicht: 

— mit internationalem Recherchenbericht 

Zur Erklarung der Zweibuchstaben- Codes und der anderen Ab- 
kurzungen wird auf die Erklarungen ("Guidance Notes on Co- 
des and Abbreviations ") am Anfang jeder regular en Ausgabe der 
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zwischen Austenit und der weiteren Phase vorliegen. Ferner wird durch Losen eines Stefan-Problems ermittelt, ob und wie sich die 
Konzentrationsverteilung (K) des beweglichen Legierungselements im austenitischen Bereich des betrachteten Volumenelements 
(9) andert und ob und um welches Ausmass (8x', 8x") sich die Phasengrenze (11, 12) dadurch verschiebt Anhand der Lage der 
Phasengrenze (11, 12) werden dann die lokalen Anteile (pi, p2, p3) der Phasen ermittelt 
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Be s chr eibung 

Rechnergesttitztes Modellierverf ahren fur das Verhalten eines 
Stahlvolumens mit einer Volumenoberf lache 

5 

Die vorliegende Erfindung betrifft ein rechnergesttitztes Mo- 
dellierverf ahren fur das Verhalten eines Stahlvolumens mit 
einer Volumenoberf lache , 

- wobei ein Rechner anhand eines momentanen Anf angszustands 
10 des Stahlvolumens und mindestens einer uber die Volumen- 

oberflache auf das Stahlvolumen einwirkenden momentanen 
Einf lussgroBe durch Losen einer Warmeleitungsgleichung und 
einer Phasenumwandlungsgleichung einen Folgezustand des 
Stahlvolumens ermittelt , 
15 - wobei die mindestens eine Einf lussgrofle fur eine Anzahl von 
Flachenelementen der Volumenoberf lache jeweils einen loka- 
len Einfluss umfasst und die lokalen Einflusse uber das je- 
weilige Flachenelement auf das Stahlvolumen einwirken, 

- wobei der Anf angszustand und der Folgezustand fur eine An- 

2 0 zahl von Volumenelementen des Stahlvolumens jeweils lokale 

Anteile von modellierten Phasen des Stahls und eine einen 
lokalen Energieinhalt des Stahls beschreibende GroBe umfas- 
sen, 

- wobei die modellierten Phasen des Stahls Austenit und eine 
25 erste weitere Phase umfassen, in die Austenit umwandelbar 

ist und die in Austenit umwandelbar ist* 

Ein derartiges Modellierverf ahren ist z. B. aus der DE-A-101 
29 565 bekannt. Bei dieser Schrift wird insbesondere erstmals 

3 0 versucht, die Fouriersche Warmeleitungsgleichung selbst und 

nicht eine unkorrekte Abwandlung dieser Warmeleitungsglei- 
chung zu losen, urn das thermodynamische Verhalten eines 
Stahlbandes korrekt zu beschreiben. Diese Schrift wird hier- 
mit durch Verweis auf sie in den Of f enbarungsgehalt der vor- 
35 liegenden Erfindung mit auf genommen . 
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In der alteren, zum Anmeldezeitpunkt der vorliegenden Erfin- 
dung nicht vorverof f entlichten deutschen Patentanmeldung 102 
51 716.9 1st ebenfalls ein derartiges Modellierverf ahren be- 
schrieben. Bei diesem Modellierverf ahren wird versucht, die 
5 Phasenumwandlung des Stahls anhand der Gibbsschen freien 

* 

Enthalpien des Stahls zu modellieren. Auch diese Schrift wird 
hiermit durch Verweis auf sie in den Of f enbarungsgehalt der 
vorliegenden Erfindung mit auf genommen . 

10 Weiterhin ist ein derartiges Modellierverf ahren z. B. aus dem 
Aufsatz "Numerische Simulation der Warmeleitung in Stahlble- 
chen - Mathematik hilft bei der Steuerung von Kuhlstrecken" 
von W. Borchers et al. r verof f entlicht ±m Unikurier der 
Friedrich-Alexander-Universitat Erlangen-Nurnberg, Heft 102, 

15 Oktober 2001/ 27. Jahrgang, bekannt . 

Schlieftlich sind noch die traditionellen Ansatze gemafi der 
Scheilschen Regel, nach Johnson-Mehl-Avrami und Brimacombe 
bekannt . 

20 

Bereits die exakte Modellierung des zeitlichen Temperaturver- 
laufs von Stahl bei der Abkuhlung, insfoesondere von Stahlban- 
dern f ist entscheidend fur die Steuerung der erf orderlichen 
Wasser- bzw . Kuhlmittelmengen einer Kuhlstrecke fur Stahl. 

2 5 Denn die beim Abkuhlen auftretende Stahlumwandlung beein- 

flusst entscheidend das thermische Verhalten des Stahls bei 
der Abkuhlung. Auch werden durch den Abkuhlprozess wesentli- 
che Materialeigenschaften des Stahls beeinflusst. Da die Ab- 
kuhlung aber nicht im thermodynamischen Gleichgewicht er- 

3 0 folgt, ist es nicht moglich, die Umwandlung einfach durch ge- 

eignete Anpassung der Warmekapazitat zu beschreiben. Daher 
ist zusatzlich eine exakte Modellierung der Phasenumwandlung 
des Stahls erf orderlich, urn die Kuhlstrecke korrekt steuern 
zu konnen. 

35 

Die traditionellen Ansatze des Standes der Technik arbeiten 
in der Praxis nicht in alien Fallen fehlerfrei. Insbesondere 
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weisen sie eine Reihe von systematischen Nachteilen auf . Zum 
Ersten muss jedes Material separat parametriert werden, In- 
terpol at ionen zwischen verschiedenen Materialien sind nicht 
Oder zumindest nur eingeschrankt moglich. Zum Zweiten werden 
5 bei den traditionellen Verfahren des Standes der Technik nur 
zwei Phasen betrachtet . Eine Erweiterung auf mebr als zwei 
Phasen ist systembedingt nicht moglich. Zum Dritten liefern 
die Verfahren des traditionellen Standes der Technik nur bei 
einer vollstandigen Umwandlung des betrachteten Metalls eine 

10 gute Ubereinstimmung zwischen Modell und Realitat. Zum Vier- 
ten liefern die Verfahren des traditionellen Standes der 
Technik keine Information uber die bei der Phasenumwandlung 
freigesetzte Umwandlungswarme . Die Kenntnis der Umwandlungs- 
warme ist fur ein korrektes Losen der Warmeleitungsgleichung 

15 aber unabdingbar. 

Die Vorgehensweisen gemaB der DE-A-101 2 9 5 65 und dem Fach- 
aufsatz "Numerische Simulation . . . " stellen dem gegenuber be- 
reits einen erheblichen Fortschritt dar, da sie zumindest die 

2 0 Warmeleitung vollig korrekt beschreiben. Die altere deutsche 

Patentanmeldung verbessert daruber hinaus die Modellierung 
der Phasenumwandlung. Insbesondere liefert sie bereits die 
bei der Phasenumwandlung auftretende Umwandlungswarme, Auch 
diese Verfahren sind aber noch verbesserungsf ahig. 

25 

Die Aufgabe der vorliegenden Erfindung besteht darin, ein Mo- 
dellierverf ahren fur ein Met all zu schaffen, das bessere Mo- 
de llergebni s s e liefert . 

3 0 Die Aufgabe wird dadurch gelost, 

— dass der Anf angszustand und der Folgezustand fur mindestens 
eines der Volumenelemente auch eine lokale Konzentrations- 
verteilung mindestens eines im Stahl beweglichen Legie- 
rungs element s umf assen, 
35 — dass im Rahmen der Umwandlungsgleichung fur das mindestens 
eine Volumenelement ermittelt wird, welche Konzentrationen 
des mindestens einen beweglichen Legierungselements beid- 



3 



WO 2005/076092 



PCT/EP2004/053709 



seits einer ersten Phasengrenze zwischen Austenit und der 
ersten weiteren Phase vorliegen, 

- dass durch Losen eines ersten Stefan-Problems ermittelt 
wird, ob und wie sich die Konzentrationsverteilung des min- 

5 destens einen beweglichen Legierungselements iin austeniti- 

schen Bereich des betrachteten Volumenelements andert und 
ob und urn welches AusmaB sich die erste Phasengrenze da- 
durch verschiebt, und 

- dass die lokalen Anteile der Phasen anhand einer durch das 
10 Ausmafi der Verschiebung der ersten Phasengrenze bestimmten 

Lage der ersten Phasengrenze ermittelt werden. 

Das in Stahl bewegliche Legierungselement ist in der Regel 
Kohlenstoff . Alternativ oder zusatzlich ist es aber auch mog- 
15 lich, dass das Legierungselement Stickstoff ist. Es kann auch 
ein anderes Legierungselement sein, das in Stahl bevorzugt an 
Zwischengitterplatzen angeordnet ist . 

Bereits dieses Modellierverf ahren ermoglicht gegenuber dem 

2 0 bisher bekannten Stand der Technik einen deutlichen Fort- 

schritt. Denn mit dieser Vorgehensweise kann. - je nach dem f 
ob die weitere Phase Ferrit oder Zementit ist - das Umwand™ 
lungsverhalten von Austenit in Ferrit bzw. von Austenit in 
Zementit und umgekehrt sehr realitatsgetreu modelliert wer- 
25 den. 

Vorzugsweise umfassen die modellierten Phasen des Stahls auch 
eine zweite weitere Phase, in die Austenit umwandelbar ist 
und die in Austenit umwandelbar ist. Denn dann ist es insbe- 

3 0 sondere moglich, 

- dass fur das betrachtete Volumenelement im Rahmen der Um- 
wandlungsgleichung auch ermittelt wird, welche Konzentrati- 
onen des mindestens einen beweglichen Legierungselements 
beidseits einer zweiten Phasengrenze zwischen Austenit und 

35 der zweiten weiteren Phase vorliegen r 

- dass durch zusatzliches Losen eines zweiten Stefan-Problems 
ermittelt wird, ob und wie sich die Konzentrationsvertei- 
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lung des mindestens einen foeweglichen Legierungselement s im 
austenitischen Bereich des betrachteten Volumenelements an- 
dert und ob und urn welches AusmaJS sich die zweite Phasen- 
grenze dadurch verschiebt, 
5 - dass die Stef an-Probleme miteinander gekoppelt sind, 

- dass den Phasengrenzen Flachenma£e zugeordnet sind, 

- dass ein Anteil des der zweiten Phasengrenze zugeordneten 
FlachenmaBes an der Summe der Flachenmafte bestimmt wird und 

- dass die lokalen Anteile der Phasen auch vom Anteil des der 
10 zweiten Phasengrenze zugeordneten Flachenmaftes an der Summe 

der Flachenmafie abhangen. 

Mit dieser Vorgehensweise ist nicht nur die Umwandlung zwi- 
schen Austenit auf der einen Seite und Ferrit Oder Zementit 
15 auf der anderen Seite sehr realitatsgetreu beschreibbar, son- 
dern insbesondere auch die Umwandlung von Austenit in Perlit 
und umgekehrt. Es ist somit das gesamte wesentliche Umwand- 
lungsverhalten von Stahl in festem Zustand korrekt beschreib- 
bar. Ferner ist mit dieser Vorgehensweise ermittelbar, ob 

2 0 Perlit gebildet wird oder nicht. 

Der Anteil des der zweiten Phasengrenze zugeordneten Flachen- 
maI5es an der Summe der Flachenmafte kann derart bestimmt wer- 
den, dass die Phasengrenzen stets nebeneinander angeordnet 
25 bleiben. Alternativ ist es aber auch moglich, diesen Anteil 
derart nachzuf uhren, dass die Phasengrenzen aufeinander zu 
streben. Anhand dieses Anteils kann dann auch abgeleitet war- 
den, ob Austenit nur in die erste weitere Phase, nur ±n die 
zweite weitere Phase oder sowohl in die erste als auch in die 

3 0 zweite weitere Phase umgewandelt wird. 



Prinzipiell ist es moglich, die Ausbildung der ersten Phasen- 
grenze dreidimensional zu ermitteln. Dies fordert aber einen 
erheblichen Rechenauf wand. Vorzugsweise wird daher die Rech- 
35 nung eindimensional vorgenommen. Dies ist einfacher moglich, 
wenn das Volumenelement ein raumf ullender aristotelescher 



5 



WO 2005/076092 



PCT/EP2004/053709 



Korper, insbesondere ein Quader, ist. Es ist also vorzugs- 
weise so, 

- dass das betrachtete Volumenelement quaderf ormig ausgebil- 
det 1st unci drei Grundabmessungen aufweist, 

5 - dass die erste Phasengrenze als Rechteck mit einer ersten 
Langsseite und einer ersten Querseite ausgebildet ist und 

- dass die erste Langsseite mit einer ersten der Grundabmes- 
sungen korrespondiert, die erste Querseite parallel zu ei- 
ner zweiten der Grundabmessungen verlauft und Verschiebun- 

10 gen der ersten Phasengrenze parallel zu der dritten der 

Grundabmessungen erfolgen . 

Bei Berucksichtigung auch der zweiten weiteren Phase ist es 
dann so, 

15 - dass die zweite Phasengrenze als Rechteck mit einer zweiten 
Langsseite und einer zweiten Querseite ausgebildet ist und 

- dass die zweite Langsseite mit der ersten der Grundabmes- 
sungen korrespondiert, die zweite Querseite parallel zur 
zweiten der Grundabmessungen verlauft und Verschiebungen 

2 0 der zweiten Phasengrenze parallel zur dritten der Grundab- 

messungen erfolgen . 

Wenn die Summe der Querseiten der Phasengrenzen ungefahr 
gleich dem 1,5- bis 3-fachen eines kritischen Lamellenabstan— 
25 des ist, bei dem eine Energiebilanz, die einerseits mit dem 

Verschieben der Phasengrenzen korrespondierende Phasenumwand— 
lungen des Stahls und andererseits mit dem Verschieben der 
Phasengrenzen korrespondierende Anderungen der Flache einer 
Grenzschicht zwischen der ersten und der zweiten weiteren 

3 0 Phase berucksichtigt , ausgeglichen ist, kann mittels des Mo- 

dells auch der Lamellenabstand des Perlits ermittelt werden. 
Der Lamellenabstand (das heiftt die Summe der Querseiten der 
Phasengrenzen) kann dabei insbesondere in etwa gleich dem 
Doppelten des kritischen Lamellenabstands sein. 

35 

Trotz der nur eindimensionalen Rechnung liefert das erfin- 
dungsgemaJie Modellierverf ahren sehr realitatsgetreue Ergeb- 
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nisse, wenn der Anteil an Austenit anhand. einer nicht linea- 
ren Funktion der Lage der Phasengrenze bzw. der Phasengrenzen 
ermittelt wird. 

< 

5 Wenn die Ermittlung, in welchen Konzentrationen das mindes- 
tens eine bewegliche Legierungselement beidseits der ersten 
Phasengrenze bzw, beidseits der ersten und beidseits der 
zweiten Phasengrenze vorliegt, anhand der Gibbsschen freien 
Enthalpien der Phasen erfolgt, ist die Ermittlung der Kon- 
10 zentrationen des mindestens einen beweglichen Legierungsele- 
ments an den Phasengrenzen besonders einfach. 

Es ist moglich, stets die gekoppelten Stef an-Probleme zu 16- 
sen und daraus abzuleiten, welche Phasen in welchem Umfang 

15 gebildet werden. Jn manchen Fallen ist dies sogar unumgang- 

lich. Manchmal ist es aber moglich, vor dem Losen der Phasen- 
umwandlungsgleichung anhand der im Anf angszustand bereits 
vorhandenen Phasen und anhand der Gibbsschen freien Enthal- 
pien vorab zu ermitteln, ob sowohl Austenit als auch die 

2 0 erste weitere Phase vorliegen bzw. ob zusatzlich zu Austenit 
und der ersten weiteren Phase auch die zweite weitere Phase 
vorliegt . 

Im Einzelf all kann die Anzahl von Volumenelementen klein 
25 sein. Im Extremfall ist sie gleich Eins. In aller Regel um- 
fasst das Stahlvolumen aber eine Vielzahl von Volumenelemen- 
ten. Wenn das Stef an -Problem bzw. die Stef an-Probleme nur fur 
einen Teil der Volumenelemente gelost wird bzw. werden und 
die lokalen Anteile der Phasen der anderen Volumenelemente 
30 anhand der lokalen Anteile der Phasen des Teils der Volumen- 
elemente ermittelt werden, kann daher die zum Modellieren des 
Verhaltens des Stahlvolumens erf orderliche Rechenleistung 
deutlich reduziert werden, ohne die Aussagekraft der Modell- 
rechnung wesentlich zu verschlechtern . Die Warmeleitungsglei- 
35 chung hingegen wird in aller Regel individuell fur jedes Vo- 
lumenelement gelost . 
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Das Modellierverf ahren kann alternativ online und in Echtzeit 
oder aber offline ausgefuhrt werden. 

Es ist z. B . moglich, 
5 - dass dem Rechner ein Erstzustand und zumindest eine ge~ 
wiinschte EndgroBe vorgegeben werden, 

- dass das obige Modellierverf ahren iterativ angewendet wird, 

- dass der Anf angszustand der ersten Iteration dem Erstzu- 
stand und der Anf angszustand jeder weiteren Iteration dem 

10 unmittelbar zuvor ermittelten Folgezustand entspricht und 

- dass anhand des nach einer letzten Iteration ermittelten 
Folgezustands eine erwartete EndgroBe ermittelt und mit der 
gewunschten EndgroBe verglichen wird. 

15 In diesem Fall kann das erf indungsgemaBe Modellierverf ahren 
alternativ online und in Echtzeit oder offline ausgefuhrt 
werden . 

Insbesondere bei einer of f line-Ausf uhrung ist es dabei mog- 

2 0 lich f dass die Einf lussgrofien der Iterationen in ihrer Ge- 

samtheit einem Einf lussgrolienverlauf entsprechen r der Rechner 
auf Grund des Vergleichs der erwarteten EndgroBe mit der ge- 
wunschten EndgroBe den Einf lussgrolienverlauf variiert und 
ausgehend vom Erstzustand das obige Modellierverf ahren erneut 
25 ausfuhrt, bis zumindest die erwartete Endgrofte mit der ge- 
wunschten Endgrofie korrespondiert. 

Es ist aber auch moglich, dass der Rechner anhand> einer aus 
dem Anf angszustand ermittelten Anf angsgrofte und einer ge- 

3 0 wunschten FolgegroBe die Einf lussgroBe ermittelt und eine Be- 

einf lussungseinrichtung derart ansteuert, dass das Stahlvolu- 
men entsprechend der ermittelten Einf lussgroBe beeinflusst 
wird. In diesem Fall muss das Modellierverf ahren online und 
in Echtzeit ausgefuhrt werden. Das Ansteuern der Beeinflus- 
35 sungseinrichtung kann dabei alternativ sofort oder in der 
darauf f olgenden Iteration erfolgen. 
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Weitere Vorteile und Einzelheiten ergeben sich aus der nach- 
folgenden Beschreibung eines Ausfuhrungsbeispiels in Verbin- 
dung mit den Zeichnungen* Dabei zeigen in Prinzipdarstellung 

5 FIG 1 eine Kuhlstrecke fur ein Stahlband, 
FIG 2 das Stahlband von FIG 1 im Detail, 

FIG 3 ein Volumenelement in perspektivischer Darstellung, 
FIG 4 ein Ablauf diagramm, 
FIG 5 ein weiteres Ablauf diagramm, 
10 FIG 6 eine Warmeleitungsgleichung und eine Phasenumwand- 

lungsgleichung, 

FIG 7 eine weitere Warmeleitungsgleichung und eine Phasen— 

umwandlungsgleichung, 
FIG 8 ein weiteres Ablauf diagramm, 
15 FIG 9 die Gibbsschen freien Enthalpien der Phasen des 

Stahls als Funktion des Kohlenstof f anteils bei einer 

ersten Temperatur, 
FIG 10 die Gibbsschen freien Enthalpien als Funktion des 

Kohlenstof f anteils des Stahls bei einer zweiten Tem- 
2 0 peratur, 

FIG 11 ein Volumenelement in perspektivischer Darstellung, 
FIG 12 das Volumenelement von FIG 11 in der Seitenansicht , 
FIG 13 ein Volumenelement analog dem von FIG 12, 
FIG 14 ein Volumenelement in perspektivischer Darstellung, 

2 5 FIG 15 das Volumenelement von FIG 14 in der Draufsicht, 

FIG 16 das Volumenelement von FIG 15 im Schnitt gemaB der 

Linie XVI - XVI von FIG 15, 
FIG 17 das Volumenelement von FIG 15 im Schnitt gemaft der 

Linie XVII - XVII von FIG 15, 

3 0 FIG 18 ein weiteres Ablauf diagramm und 

FIG 19 ein Detail von FIG 15. 

Gemaft FIG 1 soli z. B. ein Stahlband 1 gekuhlt werden, damit 
es vorgegebene Endeigenschaften annimmt . Hierzu wird das 
35 Stahlband 1 z. B. in einer Kuhlstrecke 2 mit einem Kuhlmedium 
3 - in der Regel Wasser - beauf schlagt . 
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Das Stahlband 1 weist eine Bandbreite b und eine Banddicke d 
auf • Es weist ferner eine - prinzipiell beliebige - Bandlange 
auf . Es durchlauft die Kuhlstrecke 2 mit einer Bandgeschwin- 
digkeit v. 

5 

Die Bandgeschwindigkeit v wird kontinuierlich erfasst und ei- 
nem Rechner 4 zugefuhrt r der die Kuhlstrecke 2 steuert. Der 
Rechner 4 ist somit in der Lage, in an sich bekannter Weise 
eine Wegverf olgung des Stahlbandes 1 zu realisieren. Dem 
10 Rechner 4 ist also bekannt, auf welche Bereiche des Stahlban- 
des 1 eingewirkt wird, wenn das Stahlband 1 mittels Aufbrin- 
gungseinrichtungen 5 der Kuhlstrecke 2 mit dem Kuhlmedium 3 
beaufschlagt wird. 

15 Der Rechner 4 ist mit einem Computerprogramm 6 programmiert, 

das dem Rechner 4 uber einen Datentrager 1, z . B. eine CD-ROM 
7, zugefuhrt wird. Auf dem Datentrager 7 ist das Computerpro- 
gramm 6 in (ausschlieftlich) maschinenlesbarer Form gespei- 
chert . Das Computerprogramm 6 wird vom Rechner 4 ubernommen 

2 0 und in einem Massenspeicher 8, z. B. einer Festplatte 8, des 

Rechners 4 hinterlegt . Dadurch ist der Rechner 4 in der Lage, 
bei Aufruf des Computerprogramms 6 ein Modellierverf ahren fur 
das Stahlband 1 bzw. fiir dessen einzelne Bereiche (= das 
Stahlvolumen 1) auszufuhren f das nachfolgend naher beschrie- 
25 ben wird. 

Gemaft FIG 2 wird hierzu das Stahlvolumen 1 innerhalb des 
Rechners 4 in Volumenelemente 9 zerlegt. Soweit ein Volumen- 
element 9 dabei nicht allseitig von anderen Volumenelementen 

3 0 9 umgeben ist, sind dem jeweiligen Volumenelement 9 ein oder 

zwei Flachenelemente 10 zugeordnet. Die Flachenelemente 10 
bilden in ihrer Gesamtheit eine Volumenoberf lache des Stahl- 
bandes 1 bzw. des betrachteten Bereichs des Stahlbandes 1. 

35 Der Rechner 4 realisiert auf Grund der Programmierung mit dem 
Computerprogramm 6 insbesondere ein Modell, in dem fur jedes 
Volumenelement 9 die thermischen Kopplungen mit seiner Umge- 
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bung berucksichtigt sind. Jedes Volumenelement 9 ist dabei - 
siehe FIG 3 - quaderformig ausgebildet. Es weist also drei 
Grundabmessungen A, B, C auf, die in der Regel parallel zur 
Richtung der Bandgeschwindigkeit v, zur Bandbreite b und zur 
5 Banddicke d orient iert sind. 

In einer ersten Ausgestaltung des erf indungsgemaften Model- 
lierverf ahrens wird dem Rechner 4 gemaft FIG 4 in einem 
Schritt SI ein Erstzustand Z zugef uhrt . Der Erstzustand Z urn- 
10 fasst fur jedes Volumenelement 9 des Stahlbandes 1 zunachst 

lokale Anteile pi, p2, p3 von modellierten Phasen des Stahls. 
Die Phasen konnen insbesondere Ferrit (Anteil pi) , Zementit 
(Anteil p2) und Austenit (Anteil p3) sein. 

15 Der Erstzustand Z umfasst weiterhin fiir jedes Volumenelement 
9 des Stahlbandes 1 eine Groiie H, die einen lokalen Energie- 
inhalt des Stahls des Volumenelements 9 beschreibt, Bei- 
spielsweise kann diese GroBe H die Enthalpie H des Volumen- 
elements 9 sein. Alternativ kommen auch die Temperatur oder 

2 0 die Entropie in Frage, 

SchlieBIich umfasst der Erstzustand Z fur mindestens eines 
der Volumenelemente 9 des Stahlbandes 1, vorzugsweise fur je- 
des Volumenelement 9 des Stahlbandes l r noch eine lokale Kon- 
25 zentrationsverteilung K mindestens eines im Stahl beweglichen 
Legierungselements . Das bewegliche Legierungselement kann 
insbesondere Kohlenstoff sein. Alternativ oder zusatzlich 
kommt z, B. auch Stickstoff in Frage . 

3 0 Dem Rechner 4 wird im Rahmen des Schrittes SI weiterhin zu- 

mindest eine gewunschte Endgrofie f ' * zugef uhrt, Gegebenen- 
f alls konnen dem Rechner 4 im Rahmen des Schrittes SI auch 
noch gewunschte Zwischengr 615 en zugefuhrt werden, so dass dem 
Rechner 4 gegebenenf alls sogar ein gewiinschter zeitlicher 
35 Verlauf der Grofte vorgegeben werden kann. 



11 



WO 2005/076092 



PCT/EP2004/053709 



GemaB FIG 4 wird nun in einem Schritt S2 zunachst ein Ein- 
flussgroBenverlauf festgelegt. Der Einf lussgr6I5enverlauf um- 
fasst fur eine Vielzahl von zeitlich unmittelbar aufeinander- 
folgenden Zeitpunkten jeweils eine Einflussgrofie W. Die Ein- 
5 flussgrofte W entspricht beispielsweise der auf das Stahlvolu- 
men 1 (also das Stahlband 1 bzw. dessen betrachteten Bereich) 
auf zubringenden Kuhlmittelmenge unter Beriicksichtigung ande- 
rer Einflusse auf das Stahlvolumen 1, wie beispielsweise von 
Transportrollen, der Warmekonvektion, der Warmeabstrahlung 
10 usw. . Sie umfasst fur die Vielzahl von Flachenelementen 10 

des St ahlvo lumens 1 (siehe FIG 2) jeweils einen lokalen Ein- 
fluss. Die lokalen Einflusse wirken dann uber das jeweilige 
Flachenelement 10 ~ auf das Stahlvolumen 1 ein* 

15 GemaB FIG 4 wird nunmehr in einem Schritt S3 ein Anfangszu- 
stand ZA des Stahlvolumens 1 gleich dem Erstzustand Z und 
eine Zeitbasis t auf einen Anf angszeitpunkt tO gesetzt. Der 
Anf angszeitpunkt tO entspricht dabei in der Regel dem Zeit- 
punkt, zu dem das betrachtete Stahlvolumen 1, also z. B . ein 

20 Abschnitt des Stahlbandes 1, in die Kuhlstrecke 2 eintritt. 

Sodann wird in einem Schritt S4 anhand des Einf lussgroBenver- 
laufs die zu dem durch die Zeitbasis t bestimmten Zeitpunkt 
auf zubringende Einf lussgroBe W ermittelt. In einem Schritt S5 

25 ermittelt der Rechner 4 dann anhand des momentanen Anfangszu- 
stands ZA des Stahlvolumens 1 und der uber die Volumenober- 
flache auf das Stahlvolumen 1 einwirkenden momentanen Ein- 
flussgrofie W einen Folgezustand ZF des Stahlvolumens 1* Der 
Rechner 4 lost dabei eine Warmeleitungsgleichung und eine 

3 0 Phasenumwandlungsgleichung* Der Folgezustand ZF umfasst die- 
selben Elemente K, pi, p2, p3, H wie der Anf angszustand ZA. 

Sodann wird in einem Schritt S6 die Zeitbasis t urn eine zeit- 
liche Schrittweite At inkrementiert . Als nachstes wird in 
35 einem Schritt S7 gepruft, ob die Zeitbasis t einen Endzeit- 
punkt tl erreicht hat. Der Endzeitpunkt tl entspricht dabei 
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in der Regel dem Zeitpunkt, zu dem das betrachtete Stahlvolu- 
men 1 die Kuhlstrecke 2 wieder verlasst. 

Wenn der Endzeitpunkt tl noch nicht erreicht ist, wird in ei- 
5 nem Schritt S8 der Anf angszustand ZA gleich dem unmittelbar 
zuvor ermittelten Folgezustand ZF gesetzt und sodann zum 
Schritt S4 zuruckgesprungen. 

Wenn hingegen der Endzeitpunkt tl erreicht wurde, wird die 
10 aus den Schritten S4 fois S8 bestehende Schleife verlassen. 
Wenn danach keine weiteren MaBnahmen ergriffen werden, wird 
nur ein sogenannter Prozessbeobachter realisiert. Vorzugs- 
weise aber wird, wie in FIG 4 dargestellt, in einem Schritt 
S9 anhand des nunmehr ermittelten Folgezustands ZA eine er- 
15 wartete EndgroBe f ' , z. B . die Temperatur Oder die Material- 
harte, ermittelt und mit der gewiinschten Endgr6J5e f ' * vergli- 
chen. Wenn die erwartete Endgrofte ±' mit der gewiinschten End- 
gro.Be £'* nicht korrespondiert , wird - die entsprechende Re- 
chenleistung vorausgesetzt - in einem Schritt S10 der Verlauf 

2 0 der Einf lussgro.Be W variiert . Sodann wird zum Schritt S3 zu- 

ruckgesprungen. Anderenfalls ist der erf orderliche Kuhlmit- 
telmengenverlauf ermittelt, so dass der Rechner 4 nunmehr die 
Kuhlstrecke 2 in einem Schritt Sll entsprechend ansteuern 
kann . 

25 

Der Vollstandigkeit halber sei dabei erwahnt, dass der Rech- 
ner 4, wenn ihm aufter der gewiinschten Endgro.Be f ' * auch ge- 
wiinschte ZwischengroJBen vorgegeben werden, selbstverstandlich 
den Kiihlmittelmengenverlauf erst bis zu der ersten gewiinsch- 

3 0 ten Zwischengrofie ermittelt, dann bis zur zweiten gewiinschten 

Zwischengr6.Be usw. , bis der gesamte Kiihlmittelmengenverlauf 
bis zur gewiinschten Endgro.Be f ' * ermittelt ist, 

GemaB den FIG 1 und 4 wird das obenstehend beschriebene Mo- 
35 dellierverf ahren online und in Echtzeit ausgef iihrt . Es ist 

aber selbstverstandlich ohne Weiteres auch eine Offline-Aus- 
flihrung moglich. Dies ist in FIG 1 dadurch angedeutet, dass 
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die Verbind-ung der Kuhlstrecke 2 mit dem Rechner 4 als offen- 
bare Verbindung angedeutet 1st, also nicht unmittelbar eine 
Ansteuerung der Kuhlstrecke 2 erfolgen muss. Weiterhin ist 
aus diesem Grund der Schritt Sll in FIG 4 nur gestrichelt 
5 eingezeichnet . 

Die Vorgehensweise gemaft FIG 4 erfordert einen sehr hohen Re — 
chenaufwand. Wenn die zur Ausfuhrung des Verfahrens gemaft FIG 
4 erf orderliche Rechenleistung nicht zur Verfugung steht, 
10 aber dennoch eine Online-Steuerung unter Nutzung des erfin- 
dungsgemalien Modellierverf ahrens erfolgen soil, wird online 
und in Eclat zeit ein Verfahren ausgefuhrt, das nachfolgend in 
Verbindung mit FIG 5 naher beschrieben wird. 

15 GemaB FIG 5 nimmt der Rechner 4 in einem Schritt S12 analog 
zum Schritt SI den Erstzustand Z entgegen . In einem Schritt 
S13 setzt der Rechner 4 analog zum Schritt S3 den Anfangszu— 
stand ZA gleich dem Erstzustand Z und die Zeitbasis t auf den 
Zeitpunkt tO. Sodann ermittelt der Rechner 4 in einem Schritt 

2 0 S14 eine gewunschte Folgegrofte f* Oder nimmt sie entgegen. 

In einem Schritt S15 ermittelt der Rechner 4 analog zum 
Schritt S9 anhand des Anf angszustands ZA eine Anf angsgrofie f . 
Anhand der im Schritt S15 ermittelten Anf angsgroBe f und der 

2 5 gewunschten FolgegroBe f* ermittelt der Rechner 4 dann in ei- 

nem Schritt SI 6 die Einf lussgroBe W. In einem Schritt S17, 
der mit dem Schritt Sll korrespondiert , steuert der Rechner 4 
schlieftlich die Kiihlstrecke 2 entsprechend der ermittelten 
Einf lussgroJBe W an. Das Ansteuern erfolgt dabei, wie in FIG 5 

3 0 dargestellt, vorzugsweise sofort. Gegebenenf alls kann aber 

auch ein Ansteuern im nachsten Zyklus erfolgen. 

Die nachf olgenden Schritte S18 bis S21 korrespondieren mit 
den Schritten S5 bis S8 von FIG 4. Bezuglich dieser Schritte 
35 SI 8 bis S21 wird daher von einer detaillierten Erlauterung 
abgesehen . 
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Das Verfahren gemaB FIG 5 wird vorzugsweise z . B. bei einer 
Regelung eines die Kuhlstrecke 2 durchlaufenden Stahlbandes 1 
auf einen Verlauf der gewiinschten FolgegroBe f* eingesetzt. 
In diesem Fall muss das obenstehend in Verbindung mit FIG 5 
5 beschriebene Verfahren selbstverstandlich individuell fur je — 
den einzelnen Abschnitt des Stahlbandes 1 ausgefuhrt werden . 
Auch muss fur jeden Abschnitt jeweils diejenige der Aufbrin- 
gungseinrichtungen 5 angesteuert werden, in deren Einflussbe- 
reich sich der betrachtete Abschnitt gerade befindet. Dies 
10 wird durch die eingangs erwahnte Wegverf olgung gewahrleistet . 
Das Verfahren gemafi FIG 5 wird in diesem Fall ferner parallel 
fur alle in der Kuhlstrecke 2 befindlichen Abschnitte des 
Stahlbandes 1 ausgefuhrt. 

15 Die jeweils gewunschte FolgegroBe f* kann dem Rechner 4 ex- 
plicit vorgegeben werden. Es 1st aber auch moglich, dass dem 
Rechner 4 lediglich die gewunschte Folgegrofte f* z. B. beim 
Austritt des Stahlbandes 1 aus der Kuhlstrecke 2 oder beim 
Erreichen eines Haspels, auf den das Stahlband 1 aufgehaspelt 

2 0 wird, vorgegeben wird. In diesem Fall ermittelt der Rechner 4 
selbsttatig anhand einer vorbestimmten Ermittlungsvorschrif t 
die gewiinschten FolgegroBen f* fur die einzelnen Iterationen. 

Die Regelung gemaft FIG 5 kann selbstverstandlich auch ander- 
2 5 weitig strukturiert werden, z. B. so, wie sie in der alteren 
DE 103 21 792.4 beschrieben ist . 

Wie erwahnt, werden in den Schritten S5 und S18 eine Warme- 
leitungsgleichung und eine Phasenumwandlungsgleichung gelost . 

30 

Ein Beispiel einer Warmeleitungsgleichung und einer Phasenum- 
wandlungsgleichung sind in FIG 6 dargestellt. Dort ist ein 
dreidimensionaler Ansatz fur die Warmeleitung getroffen. Bei 
Bandern kann dieser Ansatz in aller Regel aber auf den eindi- 
35 mensionalen Ansatz gemai5 FIG 7 vereinfacht werden, da der 

Warmefluss in Bandlangsrichtung und Bandquerrichtung in der 
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Regel vernachlassigbar 1st. 3/3x ist dabei die Ort sable it ung 
in Banddickenrichtung. 

Ersichtlich sind die Waxmeleitungsgleichung und die Phasenum- 
5 wandlungsgleichung miteinander gekoppelt . Zur Losung der War- 
meleitungsgleichung und der Phasenumwandlungsgleichung (= 
Schritt S5, S18) wird daher gemaft FIG 8 wie folgt vorgegan- 
gen: 

10 In einem Schritt S22 wird zunachst die lokale Temperatur T 

ermittelt . Die Ermittlung erf olgt anhand der Enthalpie H und 
der Anteile pi, p2, p3 der Phasen Ferrit, Zementit und Auste- 
nit. Auf Grund der geringen GroBe des betrachteten Volumen- 
elements 9 kann namlich angenommen werden, dass die lokale 

15 Temperatur T innerhalb des Volumenelements 9 konstant ist. 

Damit weisen auch die Phasen des Stahls innerhalb dieses Vo- 
lumenelements 9 dieselbe lokale Temperatur T auf. Die Enthal- 
pie H des Volumenelements 9 lasst sich somit als 

2 0 /? = plHl(T) + p2H2(T) + p3H3(T) ( 1 ) 

schreiben. Da ferner fur jede der Phasen die Enthalpien Hi, 
H2 bzw, H3 der jeweiligen Phase als Funktion der lokalen Tem- 
peratur T - gegebenenf alls bei Austenit unter Berucksichti- 
25 gung des Kohlenstof f anteils - eineindeutig bestimmt sind, 

lasst sich durch die Gleichung 1 die lokale Temperatur T ohne 
Weiteres ermitteln . 

* 

Fur die im Schritt S22 ermittelte lokale Temperatur T werden 
30 dann in einem Schritt S23 - getrennt fur die Phasen Ferrit, 
Zementit und Austenit - die Gibbsschen freien Enthalpien Gl 
bis G3 der einzelnen Phasen als Funktion des Anteils des be- 
weglichen Legierungselements ermittelt. Beispiele derartiger 
Verlaufe sind in den FIG 9 und 10 dargestellt. Bei dem Bei- 
35 spiel gemafl FIG 9 liegt die lokale Temperatur T dabei ober- 
halb der Umwandlungstemperatur von Austenit in Perlit, bei 
dem Beispiel gemaft FIG 10 darunter. 
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Es ist moglich, stets anzunehmen, dass alle drei Phasen Fer- 
rit, Zementit und Austenit vorliegen und erst spater - siehe 
die nachfolgenden Erlauterungen zu FIG 18 - zu entscheiden, 
welche Phasen vorliegen und welche Phasenumwandlungen erfol- 
5 gen. Gemafi FIG 8 wird aber zunachst in einem Schritt S24 die 
Anzahl der zu betrachtenden Phasen ermittelt. Dies geschieht 
wie folgt: 

Zunachst wird das System Ferrit- Austenit heraus gegriffen. 

10 Fur dieses System wird uberpruft, ob es bei dem im Volumen- 
element 9 vorhandenen Gesamtanteil des beweglichen Legie- 
rungselements rein ferritisch, rein austenitisch oder ge- 
mischt stabil ist und welche Phasenverteilung gegebenenf alls 
im stabilen Zustand vorliegt. Der Gesamtanteil des bewegli- 

15 chen Legierungselements kann anhand der Konzentrationsvertei- 
lung ohne Weiteres ermittelt werden . 

Die Ermittlung der stabilen Phase (n) und gegebenenf alls der 
Phasenverteilung geschieht dadurch, dass versucht wird, ein 

2 0 Minimum fur die gesamte Gibbssche freie Enthalpie G eines 

derartigen Systems zu ermitteln. Die konkrete Vorgehensweise 
ergibt sich dabei so, wie sie in den Seiten 16 bis 18 der Si- 
te ren deutschen Patentanmeldung 102 51 716.9 beschrieben ist, 
Diese Vorgehensweise liefert zugleich auch Konzentrationen kl 
25 und k3, in denen das bewegliche Legierungselement - typi- 
scherweise Kohlenstoff - an einer etwaigen Phasengrenze 11 
zwischen Ferrit und Austenit im ferritischen bzw. im austeni- 
tischen Bereich gegebenenf alls vorliegt. 

3 0 In analoger Weise werden sodann die Mischungssysteme Ferrit- 

Zementit (= Perlit) und Austenit-Zementit untersucht. Die Un- 
tersuchung des Mischungssystems Austenit-Zementit liefert da- 
bei zugleich auch Konzentrationen k2 und k4, in denen das be- 
wegliche Legierungselement an einer etwaigen Phasengrenze 12 
35 zwischen Zementit und Austenit im zementit is chen bzw. im 
austenitischen Bereich gegebenenf alls vorliegt. 



17 



WO 2005/076092 



PCT/EP2004/053709 



Die drei fur die drei Zweiphasen-Systeme ermittelten minima- 
len Gibbsschen freien Enthalpien G werden miteinander vergli- 
chen und das Zweiphasen-System mit der insgesamt minimalen 
Gibbsschen freien Enthalpie G herausgegrif f en . Wenn dieses 
5 Zweiphasen-System Austenit enthalt, lasst sich anhand der 
Gibbsschen freien Enthalpien ermitteln, ob der stabile Z la- 
st and eine oder zwei Phasen umfasst, welche Phase dies gege- 
benenfalls ist bzw. welche Phasen dies gegebenenf alls sind 
und welche Konzentrationen kl, k3 bzw. k2, k4 des beweglichen 
10 Legierungselements beidseits der Phasengrenze 11 von Austenit 
zu Ferrit bzw. der Phasengrenze 12 von Austenit zu Zementit 
jeweils vorliegen . 

Aus dena Anf angszustand ZA ist ferner bekannt, welche Phasen 
15 anfanglich zu welchen Anteilen pi bis p3 vorliegen. 

Wenn das Zweiphasen-System mit der minimalen Gibbsschen 
freien Enthalpie G Austenit enthalt, lasst sich also durch 
eine Zusammenschau der anfanglich vorliegenden Phasen und der 

2 0 als stabil ermittelten Phasen ermitteln, ob ein Reinsystem 

mit nur einer der Phasen Ferrit, Austenit und Zementit oder 
ein Mischsystem mit zwei oder sogar drei der Phasen Ferrit , 
Austenit und Zementit betrachtet werden muss. Auf Grund des 
Umstands, dass anhand der Gibbsschen freien Enthalpien G, Gl 
25 bis G3 auch ermittelbar ist, welche Phase stabil ist bzw. 

welche Phasen stabil sind, ist daruber hinaus auch die Rich- 
tung einer etwaigen Phasenumwandlung bereits bekannt. 

Wenn das Zweiphasen-System mit der insgesamt minimalen 

3 0 Gibbsschen freien Enthalpie G das Perlit system ist und der 

Anf angszustand ZA kein Austenit enthalt (p3 = 0) , liegt ein 
vollstandig umgewandeltes Gefuge, also ein Zweiphasen-System 
Ferrit-Zementit vor, das keiner Phasenumwandlung mehr unter- 
liegt. Wenn hingegen der Anf angszustand ZA Austenit enthalt 
35 (p3 > 0) und das Perlitsystem die minimale Gibbssche freie 
Enthalpie G aufweist, lasst sich nicht ohne Weiteres ent- 
scheiden, ob nun eine oder beide der Phasen Ferrit und Zemen- 
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tit gebildet werden. In diesem Fall wird daher vorlaufig an- 
genommen, dass beide Phasen gebildet werden, also ein Drei- 
phasensystem Ferrit-Zementit-Austenit zu betrachten ist. 

5 In einem Schritt S2 5 wird als nachstes gepriift, ob in dem be- 
trachteten Volumenelement 9 mehr als eine Phase vorliegt . 
Wenn dies nicht der Fall ist, erfolgt selbstverstandlich 
keine Phasenumwandlung. Dennoch wird in einem Schritt S2 6 fur 
dieses Volumenelement 9 eine Dif f usionsgleichung fur das be- 

10 wegliche Legierungselement angesetzt und gelost, um etwaige 

Konzentrationsschwankungen des beweglichen Legierungselement s 
innerhalb des betrachteten Volumenelements 9 auszugleichen. 
Diese Vorgehensweise ist insbesondere bei Austenit, bei dem 
der Gehalt an Kohlenstoff stark schwanken kann, von Bedeu- 

15 tung. Der Ansatz und die Losung einer derartigen Dif f usions- 
gleichung ist Fachleuten allgemein bekannt . Auf den Ansatz 
der Dif f usionsgleichung und deren Losung wird daher nachste- 
hend nicht naher eingegangen. 

2 0 Wenn in dem betrachteten Volumenelement 9 mehr als eine Phase 
vorliegt , wird als nachstes in einem Schritt S27 abgepruft, 
ob alle drei Phasen vorliegen. Wenn dies nicht der Fall ist, 
liegt ein Zweiphasen— System vor. In diesem Fall wird in einem 
Schritt S2 8 uberpriift, ob eine der beiden vorliegenden Phasen 

25 Austenit ist. Wenn dies nicht der Fall ist f liegt ein voll- 
standiges Perlit- System vor, also eine Schicht struktur aus 
Ferrit und Zementit. Diese Struktur ist im Wesentlichen sta- 
bil. In diesem Fall mussen daher keine weiteren Maftnahmen er- 
griffen werden. 

30 

Wenn hingegen eine der Phasen Austenit ist und somit die an- 
dere Phase Ferrit Oder Zementit, existiert eine Phasengrenze 
11 zwischen Austenit und Ferrit (siehe FIG 11 und 12) oder 
eine Phasengrenze 12 zwischen Austenit und Zementit (siehe 
35 FIG 11 und 13) . In diesem Fall wird als vereinf achende Maft- 

nahme angenommen, dass die Phasengrenze 11 bzw. 12 als Recht- 
eck ausgebildet ist, das eine Langsseite und eine Querseite 
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aufweist. Die Langsseite und die Querseite des Rechtecks kor- 
respondieren dabei mit den Grundabmessungen A und B des be- 
trachteten Volumenelements 9. Verschiebungen der Phasengrenze 
11, 12 erfolgen in diesem Fall, wie insbesondere aus FIG 11 
5 ersichtlich ist, parallel zur Grundabmessung C. 

Wie Fachleuten allgemein bekannt ist und auch aus FIG 12 er- 
sichtlich ist, ist in Ferrit nur sehr wenig Kohlenstof f ent- 
halten. An der Phasengrenze 11 zum austenitischen Bereich 
10 hingegen erfolgt ein Konzentrationssprung. Die Konzentratio- 
nen kl, k3, in denen das beweglich'e Legierungselement an der 
Phasengrenze 11 vorliegt, sind dabei die Konzentrationen kl, 
k3, die obenstehend im Schritt S24 fur das Ferrit-Austenit- 
System ermittelt wur den . 

15 

Wenn, wie in FIG 12 dargestellt, eine Umwandlung von Austenit 
in Ferrit erfolgen soli, muss der "Konzentrationsberg" des 
beweglichen Legierungselement s im austenitischen Bereich in 
der Nahe der Phasengrenze 11 abdif f undieren . Umgekehrt muss, 

2 0 wenn eine Umwandlung von Ferrit in Austenit erfolgen soli, 

das "Konzentrationstal" des beweglichen Legierungselement s im 
austenitischen Bereich in der Nahe der Phasengrenze 11 stan- 
dig aufgefullt werden. Die numerische oder analytische Losung 
dieser Aufgabe ist allgemein als Stefan-Problem bekannt. Sie 

2 5 erfolgt in einem Schritt S2 9. Das Stefan-Problem lasst sich 
im vorliegenden Fall wie folgt formulieren: 

K ' t - DK 'xx= 0 < 2) 

30 D ist dabei die - gegebenenf alls temperaturabhangige - Diffu- 
sionskonstante des beweglichen Legierungselement s im Auste- 
nit. K' ist die Konzentration K. Die Indizes t und x bedeuten 
die Ableitung nach der Zeit bzw. dem Ort in Verschieberich- 
tung der Phasengrenze 11. Die Verschieberichtung der Phasen- 

35 grenze 11 muss dabei nicht notwendigerweise die Banddicken- 
richtung sein. 
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Fur die Gleichung 2 sind die Anf angsbedingungen durch die 
ortliche Konzentrationsverteilung K' des beweglichen Legie- 
rungselements im Anf angszustand ZA definiert. Beim Losen des 
Stefan-Problems muss en dann nur noch die Randbedingungen be- 
5 trachtet werden, dass im f erristischen bzw. im austenitischen 
Bereich an der Phasengrenze 11 die Konzentrationen kl bzw. k3 
vorliegen, dass das bewegliche Legierungselement das betref— 
fende Volumenelement 9 nicht verlassen kann und dass sich die 
Verschiebung 8x' der Phasengrenze 11 gemaft der Stef an-Bedin- 
10 gung zu 

Qc3 - kl)dx'/dt = -DK X ( 3 ) 

ergibt, wobei St die beim Losen des Stefan-Problems benutzte 
15 zeitliche Schrittweite 1st. Diese kann dabei ein Bruchteil 

{1/2, 1/3, 1/4 der zeitlichen Schrittweite At oder aber 

gleich der zeitlichen Schrittweite At sein. 

Durch Losen des Stefan-Problems im Schritt S2 9 wird also er- 

2 0 mittelt, ob und wie sich die Konzentrationsverteilung K bzw. 

K' des beweglichen Legierungselements im austenitischen Be- 
reich des betrachteten Volumenelement s 9 andert . Zugleich 
wird dabei auch ermittelt , ob und um welches AusmaB 5x r sich 
die Phasengrenze 11 dadurch verschiebt. 

25 

Falls die weitere Phase nicht Ferrit, sondern Zementit 1st, 
ergibt sich ein prinzipiell gleich gearteter Losungsansatz 
und auch eine prinzipiell gleich geartete Losung. Der Unter- 
schied 1st nur, dass bei der Bildung von Zementit im Zementit 

3 0 Kohlenstoff stark angereichert 1st und im austenitischen Be- 

reich in der Nahe der Phasengrenze 12 zum Zementit Kohlen- 
stoff verarmt 1st. Auch dies ist Fachleuten allgemein bekannt 
und in FIG 13 dargestellt. Es andern sich also die Randbedin- 
gungen* Die Verschiebung Sx" der Phasengrenze 11 gemafi der 
35 Stef an-Bedingung ergibt sich in diesem Fall zu 

(k4- k2)8x"/dt = ~D1C' X ( 4 ) 
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K" ist dabei wieder die Konzentration K. 

Wenn hingegen im Schritt S2 8 festgestellt wurde, dass alle 
drei Phasen vorliegen bzw. vorliegen konnten, ergibt sich ein 
5 komplizierteres Problem. Denn dann muss der Ansatz gewahlt 
werden, dass eine Schichtstruktur existiert, die - siehe FIG 
14 und 15 - abwechselnd aus einer Ferrit schicht 13 und einer 
Zement it schicht 14 besteht. Die Schichtstruktur grenzt an ei- 
nen austenitischen Bereich 15 an. Die FIG 16 und 17 zeigen 
10 jeweils einen Verlauf der Konzentration K' bzw. K" des beweg- 
lichen Legierungselements in einer Ferrit schicht 13 bzw. ei- 
ner Zementitschicht 14 und im austenitischen Bereich 15 vor 
diesen Schichten 13 r 14. 

15 Im Rahmen des Modells wird dabei nach wie vor angenommen, 

dass die Phasengrenzen 11, 12 als Rechtecke ausgebildet sind, 
die jeweils eine Langsseite und eine Querseite aufweisen. Die 
Langsseiten korrespondieren weiterhin mit der ersten Grundab— 
messung A. Die Querseiten verlaufen nach wie vor parallel zur 

2 0 zweiten Grundabmessung B. Auch wird nach wie vor angenommen, 

dass Verschiebungen der Phasengrenzen 11, 12 parallel zur 
dritten Grundabmessung C erfolgen. 

Fur jede der Phasengrenzen 11 und 12 wird nunmehr in einem 
25 Schritt S30 ein Stefan -Problem angesetzt und gelost. Der 
Schritt S30 ist in FIG 18 det aillierter dargestellt. 

Das Stefan-Problem fiir die Phasengrenze 11 zwischen Ferrit 
und Austenit gehorcht gemafi FIG 18 , siehe dort Schritt S31, 

3 0 folgender GesetzmaBigkeit : 

K\-DK\=Ll (5) 

K' ist dabei die Konzentration des beweglichen Legierungsele- 
35 ments vor der Phasengrenze 11. 
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Das Stefan-Problem fur die Phasengrenze 12 zwischen Zementit 
und Austenit gehorcht folgender Gesetzmafiigkeit : 



K "t- DK "xx = L2 



(6) 



K" ist dabei - analog zu K' - die Konzentration des bewegli- 
chen Legierungs elements vor der Phasengrenze 12. 

LI und L2 sind Kopplungsterme . Sie sind — siehe FIG 14 — 
10 Funktionen eines Lamellenabstands 1 der Schichtstruktur und 
eines Anteils q der Zementitphase an der Schichtstruktur so- 
wie dessen zeitlicher Ableitung. Ausgeschrieben konnen sie z 
B. lauten: 



15 



Ll = 



1 



2 i 
q -q-1 



— 2 ^(1-^X3^ + 1) 



v 



J 



und 



(7) 



r 



L2 = 



2 1 
q -q-1 



48(1 -q)D 

* 4 #(3? -4) 

I 2 * 



(8) 



Wie aus den Gleichungen 5 bis 8 ersichtlich ist, sind die 
2 0 beiden Stefan-Probleme auf Grund der Kopplungsterme Ll, L2 

auf der rechten Seite der Gleichungen 5 und 6 miteinander ge 
koppelt . 

Auch in den Gleichungen 5 bis 8 stehen die Indizes t und x 
2 5 wieder fur die Ableitung nach der Zeit bzw. dem Ort in Ver- 
schieberichtung der Phasengrenze 11 bzw. 12. 



Fur die Verschiebungen Sx' , 8x" der Phasengrenzen 11, 12 gel- 
ten weiterhin die Stef an-Bedingungen (siehe Gleichungen 3 und 
30 4) • Es ist daher moglich, die Verschiebungen Sx' , Sx" gemaft 
den Gleichungen 3 und 4 in einem Schritt S32 gleichzusetzen . 
Dann ist der Anteil q derart zu bestimmen, dass die Verschie- 
bung 8x' der Phasengrenze 11 gemaft Gleichung 3 und die Ver- 
schiebung Sx" der Phasengrenze 12 gemaB Gleichung 4 denselben 
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Wert annehmen. Der Anteil q ist in diesem Fall also derart zu 
bestimmen, dass die Phasengrenzen 11, 12 stets nebeneinander 
angeordnet bleiben. 

5 Die Gleichheit der Verschiebungen 8x' , Sx" der Phasengrenzen 
11, 12 kann also dadurch erreicht werden, dass der Anteil q 
der Zementitphase an der Schichtstruktur des Perlit geeignet 
gewahlt wird. Man kann also - z . B . durch Einsetzen und Aus- 
rechnen (trial and error) - ermitteln, bei welchem Anteil q 
10 der Zementitphase die Verschiebungen 5x' , Sx" der Phasengren- 
zen 11, 12 ubereinstimmen . 

Anhand des so ermittelten Anteils q lasst sich dann auch 
bestimmen, ob tatsachlich Perlit gebildet wird oder ob nur 
15 eine der beiden Phasen Ferrit und Zementit gebildet wird. 

Denn nur wenn der Anteil q zwischen Null und Eins liegt, wird 
Perlit gebildet. Ist der Anteil q hingegen grower als Eins, 
wird ausschliefllich Zementit gebildet. Ist er umgekehrt klei- 
ner als Null, wird ausschlieftlich Ferrit gebildet. 

20 

Der Anteil q kann bei dieser Vorgehensweise aus numerischen 
Grunden beliebige Werte annehmen, also insbesondere auch 
kleiner als Null oder grower als Eins sein. Diese Werte sind 
aber physikalisch sinnlos. Der Anteil q wird daher in Schrit- 
25 ten S33 bis S3 6 gegebenenf alls entsprechend korrigiert. 



Die Vorgehensweise gemaJft Schritt 32 kann zu numerischen Prob 
lemen fiihren. Alternativ ist es daher moglich, unterschiedli 
che Verschiebungen 5x' , 8x" zuzulassen. In diesem Fall wird 
30 der Anteil q des Zementit an der Schichtstruktur gemaft 

q t =a — (9) 



in Schritten S32a und S32b derart nachgefuhrt, dass die Lagen 
35 der Phasengrenzen 11, 12 aufeinander zustreben. Die Schritte 
S32a und S32b werden dabei gemaJft FIG 18 anstelle des Schrit- 
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tes S32 ausgefiihrt. a ist eine geeignete Proportionalitats- 
konstante. Ihr Wert ist grower als Null. 

In diesem Fall wird eine mittlere Verschiebung 8x der Phasen- 
5 grenzen 11, 12 zu 

8x = (l-q)8x t +q8x" (10) 



def iniert . 

10 

Unabhangig davon, welche dieser beiden Vorgehensweisen 
(Schritt S32 oder Schritte S32a, S32b) ergriffen wird, wird 
aber in beiden Fallen erst durch die Bestimmung des Anteils q 
ermittelt und abgeleitet, welche Phasen in welchem Umfang ge- 
15 bildet werden. 



Wenn Perlit gebildet wird/ der Anteil q also zwischen Null 
und Eins liegt, muss noch der Lamellenabstand 1 bestimmt wer- 
den. Dies geschieht gema£ FIG 19 wie folgt: 

20 

Mit der Verschiebung Sx der Phasengrenzen 11, 12 und dem An- 
teil q ist bekannt, welche Phasenumwandlungen erfolgen und in 
welchen Anteilen sie erfolgen. Die Phasenumwandlungen liefern 
einen Beitrag 8E1 zur Energiebilanz . Der Beitrag SEl ist von 
2 5 dem Vo lumen abhangig, in dem die Phasenumwandlungen erfolgen. 
Es gilt also 

5£l=(3Al5x (ID 



3 0 (3 ist dabei eine vorab ermittelbare Proportionalitatskon— 
st ante . 

Es ergibt sich ferner eine Anderung der Flache einer Grenz- 
schicht 16 zwischen Ferrit und Zementit. Auch die Anderung 
35 dieser Flache liefert einen Beitrag 8E2 zur Energiebilanz. 
Dieser Beitrag SE2 ist proportional zur Anderung der Flache 
der Grenzschicht 16. Es gilt also 
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5£2 = 2yA5x (12) 

y 1st dabei wieder eine vorab ermittelbare Proportionalitats™ 
konstante. Der Faktor 2 ergibt sich dadurch, dass pro Phasen- 
5 grenze 11, 12 bzw. pro Schicht 13, 14 je eine Grenzschicht 16 
vorhanden 1st und zwei Phasengrenzen 11, 12 bzw. zwei Schich— 
ten 13, 14 betrachtet werden . 

Anhand der Gleichungen 11 und 12 lasst sich somit ein kriti- 
10 scher Lamellenabstand 1' ermitteln, bei dem die Energiebi- 
lanz, welche die beiden Beitrage 8E1 und 8E2 berucksichtigt , 
ausgeglichen 1st. Fur diesen kritischen Lamellenabstand 1' 

gilt 

15 l'=2y/p. (13) 

Der Lamellenabstand 1, also die Summe der Schicht dicken 11 
und 12, wird nun auf das ca. eineinhalb- bis dreifache dieses 
kritischen Lamellenabstand 1' gesetzt, z. B. auf das ca. 
2 0 zweif ache . Die Schicht dicken 11 bzw. 12 ergeben sich dann zu 

ll=(l~q)l und (14) 
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I2=ql (15) 

Die Schichtdicken 11, 12 sind proportional zu Flachenmafien 
Fl, F2, welche die Phasengrenzen 11, 12 aufweisen. Je eines 
der Flachenmafie Fl, F2 1st in FIG 14 durch Schraffieren ent- 
sprechend hervorgehoben . 



Wenn in einem der Schritte S29, S30 eine Verschiebung Sx' , 
8x" einer Phasengrenze 11, 12, ermittelt wurde, hat sich so- 
mit die Lage der Phasengrenze 11, 12 geandert . In diesem Fall 
wird in einem Schritt S37 (selbstverstandlich unter Beachtung 
35 des Vorzeichens der mittleren Verschiebung Sx) eine neue Lage 
der Phasengrenzen 11, 12 oder 11 und 12 bestimmt. In einem 
Schritt S3 8 wird dann anhand einer nicht linearen Funktion 
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der Lage der Phasengrenze 11, 12 bzw. der Phasengrenzen 11, 
12 der Anteil p3 des Austenits ermittelt. Die nicht lineare 
Funktion beriicksichtigt dabei insbesondere den Umstand, dass 
das Stefan-Problem des Schrittes S2 9 bzw. die Stef an-Probleme 
5 des Schrittes S30 eindimensional angesetzt und gelost wurden, 
in der Realitat aber eine dreidimensionale Umwandlung er- 
f olgt . 

In einem Schritt S3 9 werden sodann die Anderungen der Anteile 
10 pi/ p2 der beiden anderen Phasen Ferrit und Zementit be- 
st immt . Falls dabei die Schritte S37 bis S39 vom Schritt S29 
aus erreicht werden, ist selbstverstandlich der Anteil q Null 
oder Eins . 

15 Mit der mittleren Verschiebung 8x der Phasengrenzen 11 , 12 

und dem Anteil q des Zementit an der Schichtstruktur (0 < q < 
1) lasst sich also auf einfache Weise ermitteln, welche Ande- 
rungen sich fur die Anteile pi, p2, p3 der Phasen des Stahls 
beziiglich des betrachteten Volumenelements 9 ergeben. 

20 

' Trotz der vereinf achenden Annahmen des obigen Modellierver— 
fahrens ergibt sich ein erhebiicher Rechenaufwand. Daher wer- 
den gemaB FIG 8 vor dem Schritt S22 Schritte S4 0 und S41 und 
nach dem Schritt S39 ein Schritt S42 eingef ugt . 

25 

Im Schritt S40 werden die Volumenelemente 9 zu Gruppen zu s am- 
ine n gefasst. Beispielsweise kann eine Zusammenf assung mehre- 
rer Volumenelemente 9 erfolgen, die in Richtung der Band— 
breite b, der Banddicke d und/oder der Bandgeschwindigkeit v 

3 0 aneinander angrenzen. Auch Kombinationen hiervon sind mog- 

lich. Im Schritt S41 wird dann pro Gruppe von Volumenelemen- 
ten 9 je ein einziges der Volumenelemente 9 ausgewahlt . Nur 
fur die ausgewahlten der Volumenelemente 9 werden in den 
Schritten S22 bis S39 die Dif f erentialgleichung, das Stefan- 

35 Problem bzw. die Stef an-Probleme gelost und die Anteile pi 
p2, p3 der Phasen rechnerisch ermittelt. 
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GemsLB FIG 8 wird bei jeder Gruppe das Losungsergebnis des be- 
trachteten Volumenelements 9 im Schritt S42 fur die anderen 
Volumenelemente 9 der jeweiligen Gruppe ubernommen. Dies 
stellt die einfachste Vorgehensweise dar. Es ware aber auch 
5 denkbar, die Anteile pi, p2, p3 der Phasen von Volumenelemen- 
ten 9, deren Phasenverteilung nicht explizit berechnet worden 
ist, durch lineare oder nicht lineare Interpolation zu ermit- 
teln . 

10 Auf Grund dieses Umstands, namlich dass bei jeder Gruppe von 
Volumene lenient en 9 nur die Konzentrationsverteilung K eines 
der Volumenelemente 9 errechnet wird und fur das Losen der 
Warmeleitungsgleichung nicht die Konzentrationsverteilung K, 
sondern nur die Anteile pi, p2, p3 der vorhandenen Phasen be- 

15 notigt werden, kann gegebenenf alls auch die Konzentrations- 
verteilung K nur fur dieses eine Volumenelement 9 jeder 
Gruppe vorgegeben werden . 

Die Warmeleitungsgleichung hingegen wird in einem Schritt S4 3 
20 individuell fur jedes der Volumenelemente 9 gelost. Der 

Schritt S43 kann dabei alternativ vor oder nach dem Losen der 
Phasenumwandlungsgleichung in den Schritten S22 bis S42 aus- 
gefuhrt werden. 

2 5 Das Losen der Warmeleitungsgleichung ist ohne Weiteres mog- 

lich. Denn die Temperaturen T der einzelnen Volumenelemente 9 
sind - siehe obige Erlauterungen zum Schritt S22 — ohne Wei- 
teres ermittelbar, so dass auch deren Ortsgradient sofort er— 
mittelbar ist. Da auch die Anteile pi, p2, p3 der Phasen der 

3 0 einzelnen Volumenelemente 9 ebenfalls bekannt sind, lasst 

sich auch die Warmeleitf ahigkeit X der einzelnen Volumenele- 
mente 9 ohne Weiteres ermitteln. Da schlieBlich die Dichte p 
im Wesent lichen eine Konstante ist und die Enthalpie H der 
einzelnen Volumenelemente 9 direkt gegeben ist, ist somit 
35 auch die Warmeleitungsgleichung als Ganzes losbar. 
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Mittels der vorllegenden Erfindung wird also ein physikali- 
schen Modell erstellt, bei dem auf Basis der Gibbsschen 
freien Enthalpie G und des Dif f usionsgesetzes (Stefan-Prob- 
lem) das Umwandlungs ge s chehen hinsichtlich entstehender Pha- 
5 senanteile pi, p2, p3 und Umwandlungsgeschwindigkeit mit sehr 
hoher Genauigkeit beschreibbar ist. Es werden allgemein gul- 
tige Ergebnisse ermoglicht , welche auch die Behandlung noch 
unbekannter Materialien und Werkstoffe erlauben. Die Ansatze 
sind daruber hinaus nicht nur zur Temperaturberechnung, son— 
10 dern auch zur Berechnung von Gef iigestruktur und KorngroBe an- 
wendbar . 



Das erf indungsgemSLfre Modellierverf ahren kann an beliebiger 
Stelle eingesetzt werden, also insbesondere beispielsweise 
15 auch bei Kiihlvorgangen zwischen Walzgerusten einer WalzstraBe 
oder bei sogenanntem ferritischem Walzen. Es eignet sich auch 
zur Beschreibung der Ruckumwandlung beim Aufheizen des St ah— 
les . 
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Patentanspriiche 

1. Rechnergestutztes Modellierverf ahren fur das Verhalten ei — 
nes Stahlvolumens (1) mit einer Volumenoberflache, 
5 - wobei ein Rechner (4) anhand eines momentanen Anfangszu- 

stands (ZA) des Stahlvolumens (1) und mindestens einer liber 
die Volumenoberf lache auf das Stahlvo lumen (1) einwirkenden 
momentanen Einflussgro.Be (W) durch Losen einer Warmelei- 
tungsgleichung und einer Phasenumwandlungsgleichung einen 
10 Folgezustand (ZF) des Stahlvolumens (1) ermittelt, 

— wobei die mindestens eine Einf lussgrofte (W) fur eine Anzahl 
von Flachenelementen (10) der Volumenoberf lache jeweils ei- 
nen lokalen Einfluss umfasst und die lokalen Einflusse uber 
das jeweilige Flachenelement (10) auf das Stahlvolumen (1) 

15 einwirken, 

— wobei der Anf angszustand (ZA) und der Folgezustand (ZF) fur 
eine Anzahl von Volumenelementen (9) des Stahlvolumens (1) 
jeweils lokale Anteile (pi, p2, p3) von modellierten Phasen 
des Stahls und eine einen lokalen Energieinhalt des Stahls 

2 0 beschreibende GroBe (H) umfassen, 

— wobei die modellierten Phasen des Stahls Austenit und eine 
erste weitere Phase umfassen, in die Austenit umwandelbar 
ist und die in Austenit umwandelbar ist, 

dadurch gekennzeichnet, 
25 — dass der Anf angszustand (ZA) und der Folgezustand (ZF) fur 
mindestens eines der Volumenelemente (9) auch eine lokale 
Konzentrationsverteilung (K) mindestens eines im Stahl be- 
weglichen Legierungselement s umf as sen , 

— dass im Rahmen der Umwandlungsgleichung fur das mindestens 
30 eine Volumenelement (9) ermittelt wird, welche Konzentrati- 

onen (kl, k3; k2, k4) des mindestens einen beweglichen Le- 
gierungselement s beidseits einer ersten Phasengrenze (11, 
12) zwischen Austenit und der ersten weiteren Phase vorlie- 
gen, 

35 - dass durch Losen eines ersten Stefan— Problems ermittelt 

wird, ob und wie sich die Konzentrationsverteilung (K) des 
mindestens einen beweglichen Legierungselements im austeni- 
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tischen Bereich des betrachteten Volumenelements (9) andert 
und ob und urn welches Ausmaii (Sx, 8x' , Sx") sich die erste 
Phasengrenze (11, 12) dadurch verschiebt, und 

- dass die lokalen Anteile (pi, p2, p3) der Phasen anhand ei- 
5 ner durch das Ausmaft (Sx) der Verschiebung der ersten Pha- 
sengrenze (11, 12) bestimmten Lage der ersten Phasengrenze 
(11, 12) ermittelt werden. 

2* Modellierverf ahren nach Anspruch 1, 
10 dadurch gekennzeichnet, 

- dass die modellierten Phasen des Stahls auch eine zweite 
weitere Phase umfassen, in die Austenit umwandelbar ist und 
die in Austenit umwandelbar ist, 

- dass fur das betrachtete Volumenelement (9) im Rahmen der 
15 Umwandlungsgleichung auch ermittelt wird, welche Konzentra- 

tionen (k2 , k4; kl, k3) des mindestens einen beweglichen 
Legierungs elements beidseits einer zweiten Phasengrenze 
(12, 11) zwischen Austenit und der zweiten weiteren Phase 
vorliegen, 

2 0 - dass durch zusatzliches Losen eines zweiten Stefan -Problems 

ermittelt wird, ob und wie sich die Konzentrationsvertei— 
lung (K) des mindestens einen beweglichen Legierungsele— 
merits im austenitischen Bereich des betrachteten Volumen- 
elements (9) andert und ob und urn welches AusmaB (Sx" , Sx' ) 
25 sich die zweite Phasengrenze (12, 11) dadurch verschiebt, 

- dass die Stef an-Probleme miteinander gekoppelt sind, 

- dass den Phasengrenzen (11, 12) FlachenmaBe (Fl, F2) zuge— 
ordnet sind, 

- dass ein Anteil (q) des der zweiten Phasengrenze (12) zuge- 

3 0 ordnet en Flachenmaftes (F2) an der Summe der Flachenmafte 

(Fl, F2) bestimmt wird und 

- dass die lokalen Anteile (pi, p2, p3) auch vom Anteil (q) 
des der zweiten Phasengrenze (12) zugeordneten FlachenmaJbes 
(F2) an der Summe der Flachenmafte (Fl, F2) abhangen. 

35 

3, Modellierverf ahren nach Anspruch 2, 

dadurch gekennzeichnet, dass der An- 
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teil (q) des der zweiten Phasengrenze (12) zugeordneten Fla— 
chenmafies (F2) an der Summe der FlachenmalSe (Fl, F2) derart 
bestlmmt wird, dass die Phasengrenzen (11, 12) stets neben- 
einander angeordnet bleiben. 

5 

4. Modellierverf ahren nach Anspruch 2, 

dadurch gekennzeichnet, dass der An— 
teil (q) des der zweiten Phasengrenze (12) zugeordneten Fla- 
chenmafies (F2) an der Summe der FlachenmaBe (Fl, F2) derart 
10 nachgefiihrt wird, dass die Phasengrenzen (11, 12) aufeinander 
zu streben. 

5. Modellierverf ahren nach Anspruch 2, 3, oder 4, 
dadurch gekennzeichnet, dass anhand 

15 des Anteils (q) des der zweiten Phasengrenze (12) zugeordne- 
ten FlachenmaJ5es (F2) an der Summe der Flachenmafte (Fl r F2) 
abgeleitet wird, ob Austenit nur in die erste weitere Phase, 
nur in die zweite weitere Phase oder sowohl in die erste als 
auch in die zweite weitere Phase umgewandelt wird. 

20 

6. Modellierverf ahren nach einem der obigen Anspruche, 
dadurch gekennzeich.net, 

- dass das betrachtete Volumenelement (9) quaderf ormig ausge— 
bildet ist und drei Grundabmessungen (A, B, C) aufweist, 

2 5 — dass die erste Phasengrenze (11, 12) als Rechteck mit einer 

ersten Langsseite und einer ersten Querseite ausgebildet 
ist und 

- dass die erste Langsseite mit einer ersten der Grundabmes- 
sungen (A, B, C) korrespondiert , die erste Querseite paral— 

3 0 lei zu einer zweiten der Grundabmessungen (A, B, C) ver — 

lauft und Verschiebungen (8x f , 8x") der ersten Phasengrenze 
(11, 12) parallel zu der dritten der Grundabmessungen (A, 
B, C) erfolgen. 

35 7. Modellierverf ahren nach Anspruch 6 und einem der Anspruche 
2 bis 5, dadurch gekennzeichnet, 
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- dass die zweite Phasengrenze (12, 11) als Rechteck mit ei- 
ner zweiten Langsseite und einer zweiten Querseite ausge- 
bildet ist und 

- dass die zweite Langsseite mit der ersten der Grundabmes- 
5 sungen (A, B, C) korrespondiert , die zweite Querseite pa- 
rallel zur zweiten der Grundabmes sungen (A, B, C) verlauft 
und Verschiebungen (Sx", 8x r ) der zweiten Phasengrenze (12, 
11) parallel zur dritten der Grundabmes sungen (A, B, C) er- 
folgen. 

10 

8. Modellierverf ahren nach Anspruch 7, 

dadurch gekennzeichnet, dass die 
Summe (1) der Querseiten der Phasengrenzen (11, 12) ungefahr 
gleich dem 1,5— bis 3-fachen eines kritischen Lamellenabstan— 

15 des (1') ist, bei dem eine Energiebilanz , die einerseits mit 
dem Verschieben der Phasengrenzen (11, 12) korrespondierende 
Phasenumwandlungen des Stahls und andererseits mit dem Ver- 
schieben der Phasengrenzen (11, 12) korrespondierende Ande — 
rungen der Flache einer Grenzschicht (16) zwischen der ersten 

2 0 und der zweiten weiteren Phase berucksichtigt , ausgeglichen 
ist . 

9. Modellierverf ahren nach einem der obigen Anspruche, 
dadurch gekennzeichnet, dass das Ste- 

25 fan-Problem eindimensional angesetzt und gelost wird bzw. die 
Stef an-Probleme eindimensional angesetzt und gelost werden 
und dass der Anteil (p3) an Austenit anhand einer nicht line- 
aren Funktion der Lage der Phasengrenze (11, 12) bzw. der 
Phasengrenzen (11, 12) ermittelt wird* 

30 

10. Modellierverf ahren nach einem der obigen Anspruche, 
dadurch gekennzeichnet, dass die Er- 
mittlung, in welchen Konzentrationen (kl bis k4) das mindes- 
tens eine bewegliche Legierungselement beidseits der ersten 

35 Phasengrenze (11, 12) bzw. beidseits der ersten und beidseits 
der zweiten Phasengrenze (11, 12) vorliegt, anhand der 
Gibbsschen freien Enthalpien (Gl, G2, G3) der Phasen erfolgt. 
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11. Modellierverf ahren nach einem der obigen Anspriiche, 
dadurch gekennzeichnet, dass anhand 
der im Anfangszustand (ZA) bereits vorhandenen Phasen und an— 
5 hand der Gibbsschen freien Enthalpien (Gl, G2, G3) der Phasen 
ermittelt wird, ob sowohl Austenit als auch die erste weitere 
Phase vorliegen bzw. ob zusatzlich zu Austenit und der ersten 
weiteren Phase auch die zweite weitere Phase vorliegt. 

10 12. Modellierverf ahren nach einern der obigen Anspriiche, 
dadurch gekennzeichnet, dass das 
Stahlvolumen (1) eine Vielzahl von Volumenelementen (9) um- 
f asst f dass das Stefan— Problem bzw. die Stefan— Probleme nur 
fur einen Teil der Volumenelemente (9) gelost wird bzw. wer- 

15 den und dass die lokalen Anteile (pi/ p2, p3) der Phasen der 
anderen Volumenelemente (9) anhand der lokalen Anteile (pi, 
p2 , p3) der Phasen des Teils der Volumenelemente (9) ermit- 
telt we r den. 

2 0 13. Modellierverf ahren nach Anspruch 12, 

dadurch gekennzeichnet, dass die 
Warmeleitungsgleichung individuell fur jedes Volumenelement 
(9) gelost wird. 

25 14. Modellierverf ahren nach einem der Ansprilche 1 bis 13, 
dadurch gekennzeichnet, 

— dass dem Rechner (4) ein Erst zust and (Z) und zumindest eine 
gewunschte EndgroBe (f*) vorgegeben werden, 

— dass das Modellierverf ahren nach einem der Anspruche 1 bis 
30 10 iterativ angewendet wird, 

— dass der Anfangszustand (ZA) der ersten Iteration dem Erst- 
zustand (Z) und der Anfangszustand (ZA) jeder weiteren Ite- 
ration dem unmittelbar zuvor ermittelten Folgezustand (ZF) 
entspricht und 

35 - dass anhand des nach einer letzten Iteration ermittelten 
Folgezustands (ZF) eine erwartete Endgrofte (f 9 ) ermittelt 
und mit der gewiinschten Endgrofle (f r *) verglichen wird. 
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15. Modellierverf ahren nach Anspruch 14, 

dadurch gekennzeicb.net:, dass es on- 
line und in Echtzeit oder offline ausgefuhrt wird. 

5 

16. Modellierverf ahren nach Anspruch 14 oder 15, 
dadurch gekennzeich.net, dass die Ein- 
flussgroBen (W) der Iterationen in ihrer Gesamtheit einem 
Einf lussgroBenverlauf entsprechen, dass der Rechner (4) auf— 

10 grund des Vergleichs der erwarteten EndgroBe (f ' ) mit der ge- 
wiinschten EndgroBe (f*) den Einf lussgroBenverlauf variiert 
und ausgehend vom Erstzustand (Z) das Modellierverf ahren nach 
Anspruch 12 erneut ausfuhrt, bis zumindest die erwartete End- 
groBe (f 9 ) mit der gewunschten EndgroBe (£'*) korrespondiert . 

15 

17. Modellierverf ahren nach einem der Anspruche 1 bis 13, 
dadurch gekennzeichnet, dass es on- 
line und in Echtzeit ausgefuhrt wird, dass der Rechner (4) 
anhand einer aus dem Anf angszustand (ZA) ermittelten Anfangs- 

20 groBe (f) und einer gewunschten FolgegroBe (f*) die Einfluss— 
groBe (W) ermittelt und dass der Rechner (4) eine Beeinflus— 
sungseinrichtung (2) derart ansteuert, dass das Stahlvolumen 
(1) entsprechend der ermittelten Einf lussgroBe (W) beein— 
flusst wird. 

25 

18. Datentrager mit einem auf dem Datentrager gespeicherten 
Computerprogramm (6) zur Durchfuhrung eines Modellierverf ah- 
rens nach einem der obigen Anspruche. 

30 19. Rechner mit einem Massenspeicher (8), in dem ein Compu- 
terprogramm (6) hinterlegt ist, so dass bei Aufruf des Compu- 
terprogramms (6) von dem Rechner ein Modellierverf ahren nach 
einem der Anspruche 1 bis 17 ausfuhrbar ist. 

35 20. Beeinf lus sungseinrichtung zum Beeinflussen der Temperatur 
eines Stahlvolumens (1) , insbesondere Kuhlstrecke, 
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dadurch gekennzeichnet, dass sie von 
einem Rechner (4) nach Anspruch 19 gesteuert wird. 

21. Stahl, 

dadurch gekennzeichnet, dass er eine 
Beeinflussungseinrichtung (2) nach Anspruch 20 durchlaufen 
hat, wobei die Einf lussgroBe (W) nach Anspruch 16 oder 17 er— 
mittelt wurde. 
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